ON THE INFIMUM ATTAINED BY A REFLECTED LEVY PROCESS 
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Abstract. This paper considers a Levy-driven queue (i.e., a Levy process reflected at 0), and focuses 
on the distribution of M {t), that is, the minimal value attained in an interval of length t (where it is 
assumed that the queue is in stationarity at the beginning of the interval). The first contribution is an 
explicit characterization of this distribution, in terms of Laplace transforms, for spectrally one-sided 
Levy processes (i.e., either only positive jumps or only negative jumps). The second contribution 
concerns the asymptotics of P(Af(Tii) > u) (for different classes of functions Tu and u large); here 
we have to distinguish between heavy-tailed and light- tailed scenarios. 



1. Introduction 

The class of processes with stationary and independent increments, known as Levy processes, form 
a key object in applied probability. A substantial body of literature is devoted to Levy processes 
that are reflected at 0, sometimes also referred to as Levy-driven queues, and are regarded as a 
valuable generalization of the classical M/G/1 queues; also, the important special case of reflected 
Brownian motion is covered. 

These reflected Levy processes are defined as follows. Let X = : < e R} be a Levy process 

with (without loss of generality) zero drift: EX(1) — and X{Q) = 0. Then define the queueing 
process (or: workload process, storage process) Q = {Q{t) : t > 0} through (for c > 0) 

Q{t) := sup {X{t) - X{s) - c{t ~ s)) , 

s<t 

where it is assumed that the workload is in equilibrium (stationarity) at time 0, i.e., Q{0) =d Qe- 
We refer to this process Q as the reflection of the Levy process Y — {Y{t) : i € M} at 0, where 
Y{t) := X{t) — ct. In the sequel we normalize time such that c = 1. 

When considering the steady state Qe of the reflected process introduced above, the literature can 
be roughly divided into two categories. (A) In the first place there are results on the full distribu- 
tion of Qe, in terms of the corresponding Laplace transform. Particularly for the case of one-sided 
jumps, these transforms are fairly explicit. If X is such that it has only positive jumps and is not 
a subordinator (an increasing process), X e (which is often referred to as the spectrally positive 
case), then a generalization of the classical Pollaczek-Khintchine formula was derived [22], while 
in the case of only negative jumps, X e ,y~ (spectrally negative), Qe was seen to be exponentially 
distributed. In the Levy processes literature fl ETtl . this type of results can be found under the 
denominator fluctuation theory. We recall that there are powerful tools available for numerical in- 
version of Laplace transforms [1, 13]. (B) In the second place there are results that describe the 
asymptotics of P {Qe > u) for u large. Then one has to distinguish between results in which the 
upper tail of the Levy increments is light on the one hand, sometimes referred to as the Cramer 
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case, and results that correspond to the heavy- tailed regime on the other hand; see for instance 
ifisl] and references therein. 

In the present short communication, we consider a related problem: we analyze how long the 
process consecutively spends above a given level. More formally, we consider the distribution of 
M{t) :— infsg[o.f] Q{s), i.e., the minimum value attained by the workload process in a window 
of length t, where it is assumed that the queue is in stationarity at the beginning of the interval. 
This problem has various applications: one could for instance think of the analysis of persistent 
overload in an element of a commimication network or a node in a supply chain; see e.g. lEsI] . A 
related study on the situation of infinitely-divisible self-similar input is [2]. 

Our results correspond to both branches (A) and (B) mentioned above: in [Section 2l we present 
results on the Laplace transform of M{t), relying on known results for Levy fluctuation theory; 
we also consider the special case of Brownian motion. [Section 3l identifies the asymptotics of 
P {M{Tu) > u) for different classes of functions T„ and u large; as expected, we need to distinguish 
between heavy-tailed and light-tailed input. 
Recall that 

Y{t)^X{t)~t, Q(i) =sup(r(i)-r(s)), Af (t) = inf g(s). 

s<t se[o,t] 

For a stochastic process Z we will write: 

Z_{t) := inf Z{s), 'Z{t) := sup Z{s). 

se[oA se[o.t] 

It is well known that the process Q admits the following representation (see lEO. p. 375]): 

Q[t) = g(0) + Y{t) + max (0, -Q(0) - Y{t)) 

so that, for any u > 0, 

P (M(t) > u) = P (Q(0) + Y_{t) > u) . 

Notice that due to the independent increments property of X, the random variables Q(0) and Y_{t) 
are independent, and hence Q{0) + Y_{t) =d Qe + Y.{t)- 

2. Transforms for the spectrally one-sided case 

Let Bq denote a generic, exponentially distributed random variable with parameter q > 0, that is 
independent from the process X. In this section we evaluate the double transform, with x,q > 0, 

(1) ^{x, q) := Ee-^'*^'^') = 1 - x.Jfix, q), 

where 

^(x,g):=/ e~^"P (A/(e<,) > u) du = / e'^^P (Qe + iKe,) > u) du 

and we used the fact, that for a non-negative random variable K and any x > 0, Ee~^^ = 1 — 
X > u)du. As indicated in the introduction, we compute Jtf for Levy processes with 

one-sided jumps. 

It is worth noticing that the double transform Jiff uniquely determines the distribution of M{t). 
Determining the probability distributions from the double transform requires Laplace inversion. 
It is noted that recently, substantial progress has been made with respect to this type of inversion 
techniques. Besides the 'classical' reference ||l|], we wish to draw attention to significant recent 
progress in fl3^; the latter reference specifically addresses the multidimensional transforms, and 
also provides a fairly complete literature overview. 
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2.1. Spectrally negative case. Let tp{x) = logEe'^^^^^ denote the Laplace exponent of Y (which 
in this case is defined for every x > 0) and $((7) = sup{(? > ; ^{6) = q} be the right inverse 
of tjj. It is well known, that in this case Qe is exponentially distributed with parameter $(0) > 
0. Furthermore, the Wiener-Hopf factorisation, see [171 Theorem 6.16], yields, for any x,q > 0, 

gg2;F(e,) _ ^j-^^ 0)/k{q, x), where k{x, y) = {x — ^{y))/{^{x) — y). That is, 

Ee-i:(-«) = g x~^q) 
-tPix) - q 

Theorem 1. For a spectrally negative process X, 

^. . $(0) <^{q)+x 

^ix,q) = ^l' , , x,q>0. 

x + $(0) 

Remark 1. In the spectrally negative case, it is known that Y{eq) is exponentially distributed with 
parameter ^(g). Therefore, [Theorem II can be reformulated as 

„, / Eexp(— a;Qp) 

Eexp(-a;M(eg)) - ^ ' 



Eexp(-a;r(e<,)) 



Proof. Due to the independence of Qe and Y_{eq), and using the fact that Qe is exponentially dis- 
tributed. 



J(r{x 



,q)= f e"^" / ¥{Qe>u-z)dV{Y{eg)<z)du 

Jo J{-oo,Q] 



00 

^ — xu 



-■^mn-^)dF{Y{eq) < z)du 



e 

(-00,0] 
00 p 

e-"("+*(°»du / e*(0)MP(y(e,) < z) 
^(-00,0] 

1 -Ee'^iomM- 1 m-m 



a: + $(0) x + <i>{0) $(q) 

and the claim follows from ([1). □ 

2.2. Spectrally positive case. Let ip{x) = logEe^^'^^ be the Laplace exponent of F = —Y and 
$((7) = sup{^? > : ip{9) = q} be the right inverse of t/j. In this case, the well known Pollaczek- 
Khintchine formula gives, for x > 0, 

■ip{x) 

Furthermore, note that Y_{eq) = -Y{eq). It is known that F(e^) is exponentially distributed with 
parameter ^{q). 

Theorem 2. For a spectrally positive Levy process X, 

^ix,q) = m^^l^, X,q>0, 

Tpix) q 

where the right hand side is understood in the asymptotic sense when x = $(g), that is -0'(O)l>^(<?)'0'(<l>(g))/g^. 
Remark 2. By [17, Theorem 4.8] (understood in the asymptotic sense as well when x — ^{q)), 

l>((j) i){x) ^ q 
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Therefore. lTheorem 2l can be reformulated as 

1 exp(— a;(5e) 



lexp{—xM{eq)) = 



Eexp(-.Ty(eg))' 



Proof. First, let us assume that x ^ ^{q)- Due to the independence of Qe and Y_{eq), and using the 
fact that Y{eq) is exponentially distributed. 



oo 



je{x,q)= j e"^" / V{Y_{eq)> u- z)dV{Qe<z)du 

J [u,oo) 



/ / e-^"P(y(eq) > u- z)dMdP(Qe < ^) 
e'"^ (l - dwdP (Qe < 



[0,00) Jo 

r 

[0,00 



[0,00) 

= -/ (l-e-^^)dP(Qe < 2) 1—/ (e-*(«)^-e-"^)dP(Qe <z) 

^ 1(1 _ Ee-^'^") \ (Ee-*(«)Q= - Ee^^'^^) 

_ 1 ( ^'{'i)x\ ( $((7) X \ 



V'(x) y (x-$(g))Y q ilKx) J 

If x = ^(q), then the same computations give 

^J-(i~ ^MiM^ + V''(o)g-'&(g)V'^Wg)) ^ 

$(g) \ q ) q q 

Now the claim follows from □ 
2.3. Brownian motion. 

Theorem 3. Let X bea standard Brownian motion B = {B{t) : t e R}. Then, for each t > 0, 

P(M(t) >u)=exp(-2u) 1^2(1 + i)*(\/i) - y^exp (^^^ : 
zvhere ^'(a;) = P (A/" > x)for a standard normal random variable Af. 

Proof. Since B e Q(0) =d Qe has an exponential distribution with mean 1/2. Thus, 

P (M(t) > u) = P f Q(0) + inf (B(s) - s) > u] 

^ [ P ( inf (B(s) -s)>u-x] 2exp(-2a;)dx 
Ju Velo,t] ) 

= 2exp(-2M) / P ( sup {B{s) + s) < y ) cxp(-2y)d?/ 
Jo \se[o,i] y 

= exp(-2u)Eexp [ -2 sup (B{s) + s) ) 

V '^eio.t] y 

and the claim follows after some elementary computations (see also |0, Eqn. (1.1.3)] or 1^). □ 



ON THE INHMUM ATTAINED BY A REFLECTED LEVY PROCESS 



5 



3. ASYMPTOTICS 

In this section we consider the asjonptotics of P {M{Tu) > u) for a variety of functions r„ and u 
large. As usual, heavy-tailed and light-tailed scenarios need to be addressed separately. 

3.1. Heavy-tailed case. In this section we shall work with the following assumption about the 
Levy process X: 

Assumption 1. For a > 1, let X{1) e ^y{—a) - the class of distributions with a complementary 
distribution function that is regularly varying at oo with index ~a. Moreover, if a £ (1,2), then in 
addition 

lim —7- = p e 0, 00). 

We start with the following general proposition. 

Proposition 1. Let X bea Levy process such that EX(1) = 0. Then, for any e > 0, 

Y{u) 



lim ] 

u^oo 



1 



> e 



0. 



Proof. First note that for such X we have 
(2) 



X{t) , , Y{t) 

' a.s. or equivalently — > —1 a.s.. 



t 

Fix e > 0, then 

Y{u) 



t 



>e] = ¥[- inf Y{t) > + e 
u te[o,«] 



- inf Y{t) < -I- e] =: h{u)+h{u). 
u te[o,M] 



Now (|2]l implies that 



> -1 + £ ^ 0. 



Observe that, for any T < u. 



h{u)<V[- inf y(t) <-!-£, inf Yit) < inf Y(t) 
\ute[o,u] te[o,T] te[T,u\ 



- inf Y{t)<~l-e, inf Yit)> inf Y{t) 
u t6[o,ti] te[o,T] te[T,u\ 



< 



inf Y(t) < inf Y(t) 

te[o,T] te[T.u] 



- inf Y(t) < -l-e 
u te[T,u] 



hi{u) + l22{u). 



Again ||2j implies that hiiu) 0. As for l22{u), note that 
1 



h2{u) < 



inf X(t) < ~e 

u te[T,u] 



Now we will show, that for any 5 > and u large enough we have 

- inf X(t) < -e] < S, 

u te[T,u] J 

which proves that l22{u) and completes the proof. 
Indeed, by T can be chosen such that 

¥(yt>T : \X{t)\ <5t)>l-5. 
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Therefore, for such T and 5 < e we have 



- inf X(t) < -e 

u te\T,u\ 



- inf X(t) <~eyt>T ■.\Xit)\> 5t 

u te\T,u\ ~ I wi - 



- inf X{t) <-e,\/t>T -AXmKSt 
u te[T.u\ 



where we used the fact that the last probability equals zero. 

In the sequel we say that f{n) ^ g(n) if f{n)/g{n) 1 as n cx). 



□ 



Proposition 2. Assume that the Levy process X satisfies Assumption 1 
(i) If fin) > n, then 

P {X{n) > fin)) ^ nV (X(l) > fin)) , 



as n — > oo, n G N. 
(ii) As u ^ oo, 



a — 1 



'iXil)>u). 



Proof. Ad (i). These asymptotics can be found in, e.g., [10] for a > 2 and ii for a e (1, 2); see 
also ilj] for a recent treatment. Ad (ii). See, e.g., □ 

We now state the main result of this subsection: the exact asymptotics of P (Af (r„) > u). 



Theorem 4. Assume that the Levy process X satisfies Assumption 1 Then 

PiM{Tu)>u)r~.p{Q^>u + Tu)+TuFiX{l)>u + Tu), as u ^ oo. 

The asymptotics in [Theorem 41 can be made more explicit. Part (ii) of Proposition 2 immediately 
leads to the following corollary. 



Corollary 1. Assume that the Levy process X satisfies Assumption 1 Then 

r j^mP(X(1) > u) when T„ = o(u), 

PiM{Tu)>u)^\ |±f(yl + l)-"r„P(X(l) >7^„) when u ^ AT,,, 
[ ^T„¥iXil)>Tu) when u = o(r„), 

as u ^ oo. 

Proof of [Theorem 41 The proof consists of an upper bound and a lower bound. We use the nota- 
tion := [r„J andr+ \T„^. 

Upper bound. To prove an (asymptotically) tight upper bound for P (M(T„) > u), first we observe 
that for any e > 0, using that Q(0) —d Qe is independent of {Xit) : i > 0}, 

P (M(rj > u) < P {MiT-) >u) <¥iQ,+ XiT-) >u + T~) 

< P (Qe > (1 - e)iu + T-)) + P (X(r-) > (1 - e)iu + T')) 

+ P (Qe > eiu + T-)) P (XiT-) > sin + T')) 
=: 7r+('u) +7r+(u) +7r+(M). 



Using (i) of Proposition 2 and the strong law of large numbers for X, it is easy to show that 
TT^iu) = oin^iu)) for a fixed e. It is standard now to show that 



lim lim sup — 7-r T— 

e^O u^J P (Qe > U + Tu) 



= 1. 
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Moreover, 



1- T '^ti'^) 1 

lim iim sup „ „ , — — — — = 1 , 



due to item (i) in Proposition 2 This establishes the upper bound. 
Lower bound. As for the lower bound observe that 

P (Af (r„) >u)>¥ {MiT+) >u) >P {Qe + mi) > u,X{T+) - T+ - Y{T+) < eT+) 
> P (Qe + X(T+) >u+{l + s)T+) P {X{T+) - T+ - Y{T+) < eT+) 



By Proposition 1 (u) 1 as u — > oo. Also, 

TTf (u) > P (Qe + X{T+) - eT+/2 >u + {l + e/2)T+, X{T+) > ^eT+ /2) 

> P ({Qe >u+{l+ e/2)T+, X{T+) > -eT+ /2} U {X{T+) > u + {1 + e/2)T+}) 
= P (Qe > u + (1 + e/2)T+) P {X{T+) > -eT+ /2) + P {X(T+)) > u+{l + s/2)T+) 

- P (Qe > « + (1 + e/2)T+) P {X{T+) >u+{l + e/2)T+) 
= ■ Tr3iu)T^i{u) + (u) - 7r(r(u)7rf (u), 

where we again used that Q(0) =d Qe arid {X{t) : t > 0} are independent. By the strong law of 
large numbers, tt^ (u) 1 as it oo. Moreover, it is easy to show that tt^ {u)Tr^ (u) — o{Tr^ (u)). 
Now the lower bound follows by noting that 

hmlimmf ,^ ^ , 1, 

elQ u^oo P (Qe > U + Tu) 



and that (i) of Proposition 2 yields 

hmlimmf — , , = 1. 

elo u^oo r„p(x(i) > u + r„) 

This completes the proof. □ 



3.1.1. Stable Levy processes. Following the notation from [21], let Sa{cr, /3, fi) be a stable law with 
index a S (0,2), scale parameter a > 0, skewness parameter f3 e [—1, 1] and drift fi E M.. We 
call X an (a, /3)-stable Levy process if X is a Levy process and ^(1) has the same distribution as 

Let 



B{a,l3) := — ^ '- J 1 + 13^ tan^ i^—j sin (^— + arctan (^/3tan (^j j j 

and let X be an (a, /3)-stable Levy process with a e (1, 2) and /3 e (—1, 1]. Then, 

P(x(i)>.)^:^i^.-", 

a 

see, e.g., fii, Prop. 2.1]. NowlH 

eorem 4l can be rephrased as follows. 
Corollary 2. For an (a, (3)-stable Levy process X with a e (1,2) and /3 e (—1, 1], 

P (M(rO > u) ^ i A i)~"^(£:«r„i-" w ^ yir„, 

[ 2^31 listen i/ = o(r„). 
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3.2. Light-tailed case. In this subsection, we consider the Hght-tailed situation, also frequently 
referred to as the Cramer case. Throughout, with 0(t9) := logEexp('i9X(l)) denoting the cumulant 
function, we impose the following assumption. 

Assumption 2. Let 

/3* := sup{/3 : Ee^^^^' < oo} 
Assume that /3* > and there exists ■&* e (0, (3*) such that = i?*. 

For r > 0, define 

/(r) := sup (i9r - . 

i3>0 



Proposition 3. Under Assumption 2\ the following statements hold. 

(i) As u ^ oo, 

logP(Qe >u) d*u. 

(ii) For all u > 0, 

F{Qe>u) <e-^*''. 
(Hi) The function I obeys 

/(I) < oo and /'(I) < 
(iv) For any e > 0, 

liminf ilogP(r(M) > -eu) > -/(I) 

u— >-oo u 

Proof. For (i) and (ii), we refer to IJ, Thms. 5.1 and 5.2 of Chapter XIII] or lll6l. Proposition 5.1 and 
Remark 5.1]. For (iii), notice that /(I) = sup^>q(i? — ip{'&)) is attained for ■& s (0, d*); therefore also 
/'(I) < As for (iv), observe that 

P {Y{u) > -eu) = P (li^ e A, 



where 

A, := {f e D[Q,1] : fit) > -e,yt e [0,1]} 

and D[0, 1] is the space of cadlag functions on [0, 1]. Using sample-path large deviations results 
for Levy processes (see 1 11, Theorems 5.1 and 5.2]), we now obtain that 

liminf ilogP(i:(w) > -eu) > - inf {V^/} : /e A,nC[0,l]}, 

where ij}{f) :— /(/'(<) + 1) d<. Now observe that the path /* = is in Ag.. The stated follows by 
realizir\g that tpif) = /(I). □ 

Now we can proceed with the main result of this subsection. 



Theorem 5. Assume that the Levy process X satisfies Assumption 2 Then 



logP(M(r„) > m) u'd*-TuI{l), as u ^ oo. 

The asymptotics in lTheorem 5l can trivially be made more explicit by comparing both exponential 
decay rates. The intuition behind the following corollary is that, in large deviations language, 
the most likely path corresponding to the rare event under study first builds up from an empty 
system to level u (at time 0), and then remains at level u for the nest T„ time units; both parts of 
the path result in both contributions to the decay rate (i.e., —ui)* and — Tu/(1)). Then, depending 
on whether T„ is small or large with respect to u, one of these two contributions dominates. 
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Corollary 3. Assume that the Levy process X satisfies Assumption 1 Then 

-uir 



¥{M{Tu)>u) 



when Tu — o(u), 
-TuiA-d* + /(I)) when u - AT„, 
-T„/(l) when u = o(r„), 

as u ^ oo. 

Proof of iTheorem 5l The proof again consists of two bounds. 

Lower hound. Observe that the probability of interest is, for any e > 0, bounded from below by 

^{Q,>u + eTu) P {Y{Tu) > -eTj . 



Now the lower bound follows by combining parts (i) and (iv) of Proposition 3 and then sending 
e i 0. 

Upper hound. Observe that 

where I{t) denotes the indicator fimction l{x(t)/te(o,i)}- Thus, 
P {M{Tu) > u) < P (Qe + (X(Tu) - r„)I(T„) > u) 



l{Tu) < X 



'{Q, >u-xT^ + I{Tu))d¥ 

= I V{Q,>u~xTu + Tu)AV 
Jo 



Tu 
X{Tu) 
Tu 



X{Tu) 
Tu 



l{Tu) < X 



where the last inequality follows from part (ii) of Proposition 3 The sequence {X {u)I{u) / u} sat- 
isfies the large deviations principle on ((0, 1), B{0, 1)) with rate u and rate function /(•). Thus, 
Varadhan's lemma Theorem. 4.3.1] implies 



lim — log 

u^ca 1 ,,, 



-i>*T„(l-2;)jp 



X{Tu) 
T, 



l{Tu) < X 



inf 

xe{o,i) 



where the last equality is due to part (iii) of Proposition 3 and convexity of /(•). 



□ 
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